for k=1:100
    for i=1:100
        for j=1:100
            A(i,j) = (i*cos(k/100*pi)+j*sin(k/100*pi))/200;
            B(i,j) = (i*cos(k/100*pi)-j*sin(k/100*pi))/200;
        end
    end
    [U V X C S] = gsvd(A,B);
    imagesc(vertcat(horzcat(A,B),horzcat(U,V),horzcat(C,S)));
    colormap(flipud(gray));
    axis equal;
    M1(k) = getframe;
end

movie(M1,1,10);